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A microcanonical treatment of Hagedorn systems, i.e. finite mass hadronic resonances with an 
exponential mass spectrum controlled by the Hagedorn temperature Th, is performed. We show 
that, in the absence of any restrictions, a Hagedorn system is a perfect thermostat, i.e. it imparts 
its temperature Th to any other system in thermal contact with it. We study the thermodynamic 
effects of the lower mass cut-off in the Hagedorn mass spectrum. We show that in the presence 
of a single Hagedorn resonance the temperature of any number of Nb Boltzmarm particles differs 
only slightly from T H up to the kinematically allowed limit N B ln . For N B > N B ln however, the 
low mass cut-off leads to a decrease of the temperature as Nb grows. The properties of Hagedorn 
thermostats naturally explain a single value of hadronization temperature observed in elementary 
particle collisions at high energies and lead to some experimental predictions. 
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I. INTRODUCTION 

The statistical bootstrap model (SBM) Q, Q gave the 
first evidence that an exponentially growing hadronic 
mass spectrum <7h(to) = exp[m/T#] (m /m) a for m — > 
oo (the constants m and a will be defined later) could 
lead to new thermodynamics above the Hagedorn tem- 
perature Th- Originally, the divergence of thermody- 
namic functions at temperatures T above Th was in- 
terpreted as the existence of a limiting temperature for 
hadrons. In other words, it is impossible to build the 
hadronic thermostat above Th- A few years later an 
exponential form of the asymptotic mass spectrum was 
found in the MIT bag model ^Sj and the associated limit- 
ing temperature was interpreted as the phase transition 
temperature to the partonic degrees of freedom Q . These 
results initiated extensive studies of hadronic thermody- 
namics within the framework of the gas of bags model 
(GBM) H @. The SBM with a non-zero proper vol- 
ume of hadronic bags was solved analytically by 
the Laplace transform to the isobaric ensemble and the 
existence of phase transition from hadronic to partonic 
matter (also called the quark gluon plasma, QGP) was 
shown. Since then this techni que has been used to solve 
more sophisticated versions [y, [jjj] of the SBM and other 
statistical models [TT| . 

The major achievement of the SBM is that it naturally 
explains why the temperature of secondary hadrons cre- 
ated in hadronic collisions cannot exceed Th- However, 
this result is based on two related assumptions. First, the 
grand canonical formulation for SBM is appropriate, and 
second, the resonances of infinite mass should contribute 
to thermodynamic functions. 

Very recently, using the microcanonical formulation, 
we showed [12| that in the absence of any restrictions 
on the mass, resonances with the Hagedorn mass spec- 
trum behave as a perfect thermostat and perfect chemi- 
cal reservoir, i.e. they impart the Hagedorn temperature 
Th to particles which are in thermal contact and force 
them to be in chemical equilibrium. Similar questions 



about chemi cal p roperties of the heavy resonances were 
addressed in |13| . 

Our analysis |l2j | based on very general thermody- 
namic arguments shows that it is improper to include 
any temperature other than Th into canonical and/or 
grand canonical formulations of the statistical mechan- 
ics of any system coupled to a Hagedorn thermostat. We 
also demonstrated that the Hagedorn thermostat gener- 
ates a volume independent concentration of the particles 
in chemical equilibrium with it |12| . Thus, the entire 
framework of the SBM and GBM, which is also based 
on these two assumptions, must be revisited. In other 
words, it is necessary to return to the foundations of the 
statistical mechanics of hadrons and study the role of the 
Hagedorn mass spectrum for finite masses of hadronic 
resonances above the cut-off value m , below which the 
hadron mass spectrum is discrete. Such an analysis for 
an arbitrary value of a in g#(m) was not done in |12| . 
and, will be performed here. 

This refinement is important for understanding the 
differences and similarities between A+A and elemen- 
tary particle collisions at high energies. There are two 
temperatures measured in A+A collisions that are very 
close to the transition temperature Tx r from hadron gas 
to QGP calculated from the lattice quantum chromody- 
namics |l lj at vanishing baryonic density. The first is 
the chemical freeze-out temperature at vanishing bary- 
onic density Tchem ~ 175 ± 10 MeV of the most abun- 
dant hadrons (pions, kaons, nucleons etc) extracted from 
particle multiplicities at highest SPS 15] and all RHIC 
[16| energies. Within the error bars Tr r ~ Tchem is 
also very close to the kinetic freeze-out temperature Ticin 
(i.e. hadronization temperature) found from the trans- 
verse mass spectra of heavy, weakly interacting hadrons 
such as fi hyperons, J/i/j and tp' mesons at the highest 
SPS energyJFljra, and ^-hyperons HJ E2 and 
meson [H ETISll at s/s = 130 A-GeV and y/s = 200 
A-GeV energies of RHIC. The existence of the deconfine- 
ment transition naturally explains the same value for all 
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these temperatures. 

Can the same logic be applied to the collisions of el- 
ementary particles, where the formation of a deconfined 
quark gluon matter is rather problematic? The fact that 
the hadronization temperature psl ] and inverse slopes of 
the transverse mass spectra of various hadrons [24| found 
in elementary particle collisions at high energies are sim- 
ilar to those ones for A+A collisions is tantalizing. 

In the present paper we show that the results of A+A 
and elementary particle collisions can be understood and 
explained on the same footing by the Hagedorn thermo- 
stat concept. For this purpose we study the properties of 
microcanonical equilibrium of Boltzmann particles which 
are in contact with Hagedorn thermostats and elucidate 
the effect of the mass cut-off m of the Hagedorn spec- 
trum on the temperature of a system at given energy. 
The microcanonical formulation of the Hagedorn ther- 
mostat model (HTM) is given in Section II. Section III 
is devoted to the analysis of the most probable state of a 
single heavy Hagedorn thermostat of mass m > m a and 
Nb Boltzmann particles. The last section contains our 
conclusions and possible experimental consequences. 

II. HAGEDORN THERMOSTAT MODEL 

Let us consider the microcanonical ensemble of Nb 
Boltzmann point-like particles of mass tub and de- 
generacy gB, and Njj hadronic point-like resonances 
of mass rriH with a mass spectrum <?#(to#) = 
exp[mjj/Tjj](m /mjj) a for mjj > m which obeys the 
inequalities m 3> Th and m >mj. A recent analysis 
|2fij suggests that the Hagedorn mass spectrum can be 
established for m Q < 2 GeV. 

In the SBM [H and the MIT bag model it was 
found that for m# — > oo the parameter a < 3. For finite 
resonance masses the value of a is unknown, so it will be 
considered as a fixed parameter. 

The microcanonical partition of the system, with vol- 
ume V, total energy U and zero total momentum, can be 
written as follows 
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where the quantity 
and e(M,P) 



e(m H ,Qi) (ef = e(m B ,Pj) 



\JM 2 + P 2 ) denotes the energy of the 

Hagedorn (Boltzmann) particle with the 3-momentum Qi 
(pj). In order to simplify the presentation of our idea, 
Eq. Q accounts for energy conservation only and ne- 
glects momentum conservation. 

The microcanonical partition (JTJ can be evaluated by 
the Laplace transform in total energy U 27]. Then the 
momentum integrals in Q are factorized and can be per- 
formed analytically. The inverse Laplace transform in 



the conjugate variable A can be done analytically for the 
nonrelativistic and ultrarelativistic approximations of the 
one-particle momentum distribution function 
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where the auxiliary integral is denoted as 
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Since the formal steps of further evaluation are similar 
for both cases, we discuss in detail the nonrelativistic 
limit only, and later present the results for the other case. 
The nonrelativistic approximation (MRe(X) » 1) for 
Eq. JTJ is as follows 
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where Ekin — U — muNn — msNs is the kinetic energy 
of the system. 

As shown below, the most realistic case corresponds to 
the nonrelativistic treatment of the Hagedorn resonances 
because the resulting temperature is much smaller than 
their masses. Therefore, it is sufficient to consider the 
ultrarelativistic limit for the Boltzmann particles only. 
In this case (MRe(X) « 1) the equation can be 
approximated as 
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where the kinetic energy does not include the rest energy 
of the Boltzmann particles, i.e. Eki n = U — uihNh- 

Within our assumptions the above results are general 
and can be used for any number of particles, provided 
Nh + Nb > 2. It is instructive to consider first the 
simplest case Nh = 1. This oversimplified model, in 
which a Hagedorn thermostat is always present, allows 
us to study the problem rigorously. For Nh = 1 wc treat 
the mass of Hagedorn thermostat ran as a free parameter 
and determine the value which maximizes the entropy of 
the system. The solution m* H > of 
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provides the maximum of the system's entropy, if for 
77i0 = m* H the second derivative is negative 
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If the inequality Q is satisfied, then the extremum condi- 
tion © defines the temperature of the system of (7Vb + 1) 
nonrelativistic particles 
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Thus, as to^ — > oo it follows that T*(m* H ) — ► T#, while 
for finite 3> Tjj an d a > | (a < |) the temperature 
of the system is slightly larger (smaller) than the Hage- 
dorn temperature, i.e. T* > T H (T* < T H )- Formally, 
the temperature of the system in equation JSJ may dif- 
fer essentially from Th for a light thermostat, i.e. for 
m* H < Th- However, it is assumed that the Hagedorn 
mass spectrum exists above the cut-off mass m 3> Th, 
thus m* > T H - 



III. THE ROLE OF THE MASS CUT-OFF 

Now we study the effect of the mass cut-off of the Hage- 
dorn spectrum on the inequality J7J) in more detail. For 
a < | the condition J7J) is satisfied. For a > | the in- 
equality (JZJ) is equivalent to 
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which means that a Hagedorn thermostat should be mas- 
sive compared to the kinetic energy of the system. 

A more careful analysis shows that for a negative value 
of the determinant D nr (N = Nb — |a) 

D nr = (U -m B N B -lT H iV) - 

4 (a - |) T ff (U-m B N B ) < 0, (10) 

equation © has two complex solutions, while for 
D nr — there exists a single real solution of ©. Solv- 
ing (JTOJ for (U - m B N B ), shows that for N > |a — 1, 
i.e. for N B > |, a — 1 the inequality (|1U[) does not hold 
and D nr > 0. Therefore, in what follows we will assume 
that N B > |a — 1 and only analyze the case D nr > 0. 



For this case equation © has two real solutions 
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For a < | only solution is positive and corresponds 
to a maximum of the microcanonical partition fl nr . 

For a > | both solutions of I© are positive, but only 
m H is a maximum. From the two limiting cases: 
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for E kin *0, (13) 



and the fact that m H obey the inequalities 
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Using Eq. (|11|) for m H , it is clear that for any value of 
a the constraint m H > m is equivalent to the inequality 
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Thus, at fixed energy U for all Nb < N B m at m* H 
there is a local maximum of the microcanonical par- 
tition Q nr with the temperature T = T*(m H ). For 
Nb > N B ™ the maximum of the partition Q, nr can- 
not be reached due to the cut-off constraint and, con- 
sequently, the most probable state corresponds to m B = 
to with T < T*(m ) from Eq. ©. In other words, for 
Nb > N B ln the amount of energy U is insufficient for the 
mass of the Hagedorn thermostat to be above the cut-off 
to and simultaneously maintain the temperature of the 
Boltzmann particles according to Eq. ©. By assump- 
tion there is a single Hagedorn thermostat in the system, 
therefore, as Nb grows the temperature of the system 
decreases from T*(m ) value. Thus, the equality l(TB|) 
defines the kinematical limit for reaching the maximum 
of the microcanonical partition. 

To prove that the maximum of the microcanonical par- 
tition at niH = m H is global it is sufficient to show that 
the constraint m H > to q is not consistent with the condi- 
tion m H > to . For a < | the maximum is global because 
for < van < m H (mH > m H ) the partition fl nr (NH = 
1, win) monotonically increases (decreases) with run ■ For 
a > | it is clear that the maximum at tyih — r m H is local, 
if the state with mass m B — m is more probable, i.e. 
&>nr(N H = 1,to ) > fl nr (N H = l,m H ). Due to JTJJ this 
can occur, if m H > m . Substituting Eq. (|ll|l into the 
last inequality, shows that this inequality reduces to the 
condition Nb > Ng m . This contradicts the constraint 
m H > to in the form of Eq. 115|) . Thus, the maximum 
of the microcanonical partition is global. 

To complete our consideration of the nonrelativistic 
case let us express the partition (0J in terms of the tem- 
perature JSJ). Applying the Stirling approximation to 
the factorial (|(A B + 1) - 1)! for N% in > N B > 1 
and reversing the integral representations J2Jl and 10 for 
A = l/T*(m H ), one finds 
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This is just the grand canonical partition of (N B + 1) 
Boltzmann particles with temperature T*(m H ). If Ng > 
N% m > 1, then T*(m+) in (JTSJl should be replaced by 

rp I at \ _ 2(U—mBNB—m ) 
1 (N B ) = 3 (JVb + 1) • 

Fig. 1 shows that for a > | the system's tempera- 
ture T = T*{m H ) as a function of N B remains almost 
constant for Nb < Ng, reaches a maximum at Ng ln 
and rapidly decreases like T = T (N B ) for N B > N§ n . 
For a < | the temperature has a plateau T = T*(m H ) 
for N B < N B ln , and rapidly decreases for N B > N B m 
according to T (N B ), 

The same results are valid for the ultrarelativistic 
treatment of Boltzmann particles. Comparing the non- 
relativistic and ultrarelativistics expressions for the mi- 
crocanonical partition, i.e. equations l@J and JJjJ, respec- 
tively, one finds that the derivation of the ultrarelativis- 
tic limit requires only the substitution Nb — > 2Nb and 
itlb/Th — > in equations JJ3~EEJl- Note that this substi- 
tution does not alter the expression for the temperature 
of the system, i.e. the right hand side of JSJ. 

Finally, we show that for a heavy Hagedorn thermo- 
stat (m H » m ) these results remain valid for a single 
Hagedorn thermostat split into Njj pieces of the same 
mass. Substituting m# — > tuhNh in the nonrelativistic 
expressions @ and minimizing it with respect to ran , the 
temperature of the system in the form of equation (jSJ) 
is T* (mjjNu), where the mass of Nh Hagedorn ther- 
mostats m* H is related to the solution m H of equation 
as m* H = rrijj/NH- Since the original single ther- 
mostat of mass m H was assumed to be heavy, it follows 
T*(m^N H ) = T*{m+) — > T H . A more careful study (see 
also [13) using an exact expression for the microcanon- 
ical partition of Njj Hagedorn thermostats of the same 
mass vriH gives the same result, if mn 3> m . A gen- 
eralization of these statements to the case of Njj heavy 
Hagedorn thermostats of different masses also leads to 
the same result. Thus, splitting a single heavy Hagedorn 
thermostat into an arbitrary number of heavy resonances 
(heavier than m ) does not change the temperature of the 
system. 

IV. CONCLUSIONS 

In the present work we generalized the SBM results 
[2^| to systems of finite energy by showing explicitly that 
even a single resonance with the Hagedorn mass spec- 
trum degeneracy, i.e. a Hagedorn thermostat, keeps an 
almost constant temperature close to Th for any number 
of Boltzmann particles 3 < Nb < iVjj™. For the high 
energy limit U S> m this means that a single Hagedorn 
resonance defines the temperature of the system to be 
only slightly different from Th until the energy of the 
Hagedorn thermostat is almost negligible compared to 
U. In contrast to the grand canonical formulation of the 
original SBM [2(j , in the presence of a Hagedorn thermo- 
stat the temperature Th can be reached at any energy 
density. 

The thermostatic nature of a Hagedorn system obvi- 
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FIG. 1: A typical behavior of the system's temperature as 
the function of the number of Boltzmann particles Nb for 
a — 3 and a — for the same value of the total energy U = 
30 ms. Due to the thermostatic properties of a Hagedorn 
resonance the system's temperature is nearly constant up to 
the kinematically allowed value N B ' n given by 1151 . 



ously explains the ubiquity of both the inverse slopes of 
measured transverse mass spectra 1 12 11 and hadronization 
temperature found in numerical simulations of hadrons 
created in elementary particle collisions at high energies 
|23l 123, • By a direct evaluation of the microcanonical 
partition we showed that in the presence of a single Hage- 
dorn thermostat the energy spectra of particles become 
exponential with no additional assumptions, e.g. phase 
space dominance [30| or string tension fluctuations |3l| . 
Also the limiting temperature found in the URQMD cal- 
culations made in a finite box |3^ | can be explained by 
the effect of the Hagedorn thermostat. We expect that, if 
the string parametrization of the URQMD in a box [32] 
was done microcanonically instead of grand canonically, 
then the same behavior would be found. 

The Hagedorn thermostat model generalizes the sta- 
tistical hadronization model which successfully describes 
the particle m ultip licities in nucleus-nucleus and elemen- 
tary collisions |U l28l |29| . The statistical hadronization 
model accounts for the decay of heavy resonances (clus- 
ters in terms of Refs. [23l [28l l29j l only and does not 
consider the additional particles, e.g. light hadrons, free 
quarks and gluons, or other heavy resonances. As we 
showed, the splitting of a single heavy Hagedorn reso- 
nance into several does not change the temperature of 
the system. This finding justifies the main assumption 
of the canonical formulation of the statistical hadroniza- 
tion model 28] that smaller clusters may be reduced to a 
single large cluster. Thus, recalling the MIT Bag model 
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interpretation of the Hagedorn mass spectrum |{J Q , we 
conclude that quark gluon matter confined in heavy res- 
onances (hadronic bags) controls the temperature of sur- 
rounding particles close to Th, and, therefore, this tem- 
perature can be considered as a coexistence temperature 
for confined color matter and hadrons. Moreover, as we 
showed, the emergence of a coexistence temperature does 
not require the actual deconfinement of the color degrees 
of freedom, which, in terms of the GBM ?8|, is equiva- 
lent to the formation of the infinitely large and infinitely 
heavy hadronic bag. 

Within the framework of the Hagedorn thermostat 
model we found that even for a single Hagedorn thermo- 
stat and a > § the system's temperature T = T*(m H ) 
as a function of Nb remains almost constant for Nb < 
Ng m , reaches a maximum at Ng and rapidly decreases 
for Nb > Ng n (see Fig. 1). For a < § the temperature 
has a plateau T = T*(m+) for N B < N% in , and rapidly 
decreases for Nb > Ng n . If such characteristic behavior 
of the hadronization temperature or the hadronic inverse 
slopes can be measured as a function of event multiplic- 
ity, it may be possible to experimentally determine the 
value of a. For quantitative predictions it is necessary to 



include more hadronic species into the model, but this 
will not change our result. 

If we apply the HTM to elementary particle collisions 
at high energy, then, as shown above, the temperature 
of created particles will be defined by the most probable 
mass of the Hagedorn thermostat. If the most probable 
resonance mass grows with the energy of collision, then 
the hadronization temperature should decrease (increase) 
to Th for a > | (a < §). Such a decrease is observed 
in reactions of elementary particles at high energies, see 
Table 1 in Ref. H}. 

In order to apply these results in a more physical fash- 
ion to the quark gluon plasma formation in relativistic 
nucleus-nucleus collisions (where the excluded volume ef- 
fects are known to be important 0, H, El for all 
hadrons) the eigen volumes of all particles should be 
incorporated into the model. For pions this should be 
done in relativistic fashion |34[ • Also the effect of finite 
width of Hagedorn resonances may be important |lCj| and 
should be studied. 
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